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I. INTRODUCTION 



The possibility that baryonic matter may undergo a transition to a new hadronic phase 
consisting of a condensate of negatively charged kaons {K~) in high density matter has 
been extensively discussed . The kaon condensation may be formulated on the basis 
of SU(3)lX SU(3)r current algebra and PCAC in a model-independent way 0,0] . Within 
this framework, the s-wave kaon-condensed state \6) is generated by a chiral rotation of 
the normal state |normal) as \9) = i/x|normal) , where Uk = exp{i/j,KtQ) exp{i9Q^) is a 
unitary operator. Here fix is the kaon chemical potential, 6 the chiral angle, which is the 
order parameter of the condensate, and Q (Ql) is the electromagnetic charge operator 
(the axial charge operator). The classical K~ field is then written as (K^) = {9\K^\9) = 
{f /^/2) sm9 exTp{—ijjKt) with /(=93 MeV) the meson decay constant. 

Kaon condensation is a form of Bose-Einstein condensation (BEC): the lowest exci- 
tation energy of the K~ decreases in a nuclear medium with an increase in the baryon 
number density due to an attractive s-wave kaon-nucleon [KN) interaction. When the 
lowest excitation energy becomes equal to the charge chemical potential, the distribution 
function for the K~ diverges and a macroscopic number of kaons appear. The relevant 
s-wave KN interaction comes mainly from the scalar interaction (the KN sigma term, 
Skn) and the vector interaction (the Tomozawa- Weinberg term) 0-0. 

The baryon number density for the onset of kaon condensation in stable neutron-star 
matter has been predicted to be 3—4 times the nuclear matter density uq (=0.16 fm~^), 
depending on the magnitude of Skn- Motivated by studies of kaon condensation, the 
in-medium KN interaction has recently been elaborated both theoretically and experi- 
mentally, through KN scattering 0-^], kaonic atoms and kaon-antikaon production 



via relativistic nucleus-nucleus collisions [|TlHT5|] . The experimental results suggest a sub- 



stantial decrease in the antikaon effective mass, which in turn favors the possible existence 
of a kaon condensate in the core of a neutron star. 

The appearance of the kaon-condensed phase in a neutron star would lead to soften- 
ing of the equation of state (EOS) [^,12,16| and also accelerate cooling of neutron stars 



|T^-^. The softening of the EOS associated with kaon condensation has been used to 
construct a model in which the collapse of a hot protoneutron star to a black hole is 
delayed [pT|-p^ initially, there is no kaon condensation in protoneutron stars, since the 



neutrino degeneracy and thermal effects raise the threshold density for kaon condensation. 
A kaon condensate is formed during the evolution into a compact neutron star through 
deleptonization and cooling. The maximum mass of the kaon-condensed neutron star is 



estimated to be ~1.5 Mq with Mq the solar mass p|, p!2| , p!6| , which is smaller than the 
maximum mass of ordinary neutron stars (~ 2.0Mq) |^ because of the significant soft- 
ening of the EOS due to the appearance of the condensates. If the mass of a protoneutron 
star after deleptonization or cooling exceeds the maximum mass of the kaon-condensed 
neutron star, the star eventually collapses to a black hole. Such a scenario for the forma- 
tion of a low-mass black hole (1.5—2.0 Mq) has been presented by Brown and Bethe 
to interpret the absence of a pulsar signature in SN 1987A Following their scenario. 



^Similar delayed collapse mechanisms of a hot protoneutron star to a more compact neutron 
star or to a black hole, attributed to other hadronic phase transitions, have been considered by 



some authors [p.24- 
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Baumgarte et al. made a numerical simulation for the delayed collapse of a protoneutron 
star by the use of an EOS including the phase transition to the kaon-condensed phase 

In these studies, the equilibrated EOS with kaon condensation was used. The kaons 
are created through weak interactions that change strangeness during the appearance 
and the growth of the condensate, and the time scale for the weak reactions are much 
larger than those for the strong and electromagnetic reactions which are responsible for 
the thermal equilibration of the system. In addition, there are at least three time scales 
which characterize the dynamical evolution of a newly born neutron star: the time scale 
for gravitational collapse (of order of a millisecond) and those for deleptonization and 
initial cooling (of order of a second and ten seconds, respectively). If the weak reactions 
which are responsible for the formation of a condensate proceed on a longer time scale than 
these time scales, they may control the dynamical evolution of a neutron star just after 
a supernova explosion. In this case, the nonequilibrium processes in kaon condensation 
which are brought about by the weak reactions must be considered. 

The time evolution of the system is determined by the Hamiltonian, 

H -^^strong ~l~ H-fj,=g^^ , (1) 

with the strong (weak) Hamiltonian -/^strong (-^weak)- In our case the state \6{t),t) is 
specified by the parameters, 9, fin (f^p) the neutron chemical potential (the proton chemical 
potential), fix, and /ig the electron chemical potential, for given temperature T and 
baryon number density ns- Alternatively we may choose another set of operators, 6, the 
number densities fij [i = p, n, K, e) to specify the state. If the weak Hamiltonian i/weak is 
negligibly small, then the time scale of the weak interaction is very large (r^eak ^ ^strong)- 
In such a case, it is obvious that all the number densities are time-independent, rij=0, 

which in turn implies 9=0, or they can commute with -f^strong, [^^i, -f^strong]=0, due to 
baryon number, strangeness, charge and lepton number conservations. Thus we can find 
the simultaneous eigenstate of -^strong and hi. Hence, if we apply the adiabatic (Born- 
Oppenheimer) approximation, then the time-dependence of nj(t) = {9{t),t\hi\9(t),t) is 
determined by the weak interaction Hamiltonian -ffweak, which leads to the rate equations. 
With ni{t) taken to be a function of time, it is possible to fix the time and construct an 
instantaneous eigenstate of -^strong and 9, by considering the extremum conditions for the 
thermodynamic potential Q given only by the strong interaction Hamiltonian -f^strong? 

\9{t),t) ~ \9{t),nK{t),nn{t),nj,{t),n,{t)). (2) 

Solving the rate equations then allows the temporal change of Ui to be determined, while 
the time dependence of the order parameter 9 is determined by use of the thermodynamic 
potential Q. 

The nonequilibrium processes in the kaon-condensed state \9{t),t) have been previously 
investigated in [^, with account taken of both the kaon-induced Urea process and the 



modified Urea process. The former process (abbreviated to KU-F and KU-B, where F and 
B stand for forward and backward reactions, respectively) is represented by the reactions 

n{p) n{p) + + Ue, (3a) 
n{p) + e~ ^ ■n{p) + Ue , (3b) 
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and provides the most efficient cooling mechanism for the star via neutrino and antineu- 
trino emissions [0,^. In terms of the condensate {K~), by which the system is supphed 
with energy and the reactions become kinematically possible, the KU process can be 
written symbohcally as n{p) + {K~) n{p) + e~ + z/e, n{p) + e" ^ n(p) + {K~) + z/g. 
The modified Urea process (MU-F and MU-B) represented by the reactions 

+ n — > n + p + e~ + z^e, (4a) 
n + p + e'^n + n + Ve , (4b) 

is a standard cooling process for a normal neutron star.| Because the transition matrix 
element for KU is proportional to sin 6, this process is operative only in the presence of 
a condensate. Thus, spontaneous creation of a condensate from the noncondensed state 



cannot occur solely by way of the KU and MU reactions, and, for this reason in Ref. |3T 
a small seed of the condensate {6 ^ 0) was put in by hand. 

In the noncondensed state, kaons are produced thermally by the weak reactions 

n + n ^ n + p + K~ , (5a) 
n + p + K~ ^ n + n , (5b) 

where a spectator neutron must take part in the reactions such that the kinematical 
condition in a degenerate Fermi system is satisfied. The thermal kaons are then converted 
to a condensate. Thus we can discuss the onset and the growth of a kaon condensate 
consistently without any ad hoc seed (see Sec.fT|), in which case the process (^ is primarily 
responsible for the onset of condensation. We will refer to the reaction (^) as the thermal 
kaon process, and abbreviate them to KT-F and KT-B, respectively. 

In a previous paper |^ (Paper I), we have obtained the reaction rate for KT, and 
discussed the effects of thermal kaons on the KT reaction rate in both the noncondensed 
and condensed states. We have also compared the KT reaction rate with the rates for 
the KU and MU reactions and obtained the following results: (i) In the noncondensed 
state, where the system is far from chemical equilibrium, hard thermal kaons with large 
momenta make the major contribution to the KT reaction rate, whereas in the condensed 
state, the soft thermal kaons, which appear as a Goldstone mode from the spontaneous 
breaking of V-spin symmetry, contribute significantly to the KT reaction rate, (ii) The 
KT reaction rate is larger than the rates for the KU and MU reactions over the relevant 
temperatures and baryon number densities, (iii) The KT process is dominant in the kaon- 
condensed state in chemical equilibrium as well as in the noncondensed state, and may 
determine the evolution of the system. 

Based on the above results for the KT reactions, we in this paper, investigate the 
nonequilibrium weak processes in kaon condensation by taking into account the KT pro- 
cess as well as the KU and MU processes, and clarify the effects of thermal kaons on the 
kinetics of condensation. Assuming that the system is in thermal equilibrium, we solve 
the rate equations given by the relevant weak reactions (0) — (H), and discuss the tempo- 
ral behavior of the number densities of the chemical species, the order parameter of the 



^ Throughout this paper, we do not take into account the direct Urea process in the nor- 
mal phase [34,p5l or in the kaon-condensed phase [16,pn], and also omit other weak reactions 



associated with hyperons |36|, taking the same standpoint as in Paper I. 
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condensate, and other physical quantities. In general, the energy, produced in the course 
of the nonequilibrium process, is dissipated into the surroundings, resulting in a rise in 
the temperature of the system. For simplicity, however, the temperature is kept constant 
throughout this paper. We find two characteristic time scales for the onset of a condensate 
and its subsequent growth. The magnitudes of these time scales are compared with other 
time scales characterizing the collapse of a newly-born neutron star, and possible effects 
on stellar collapse are discussed. 

The paper is organized as follows: The description of the nonequilibrium state based 
on the thermodynamic potential is addressed in Sec. |T|, and the formulation for obtaining 
the rate equations is given in Sec. ^TT|. The numerical results are then presented in Sec. 0, 
and summary and concluding remarks are given in Sec. 0. In Appendix A, an asymptotic 
behavior of the system near chemical equilibrium is discussed. Specifically, an expression 
for the relaxation time near chemical equilibrium is derived analytically. 



II. DESCRIPTION OF THE NONEQUILIBRIUM STATE 
A. Thermodynamic potential 

Here we define the physical conditions for a nonequilibrium kaon-condensed state, 
|6'()f:), The system is described by a thermodynamic potential Q. The microscopic 
quantities, such as the excitation energies of the baryons and the thermal and condensed 
kaons, are determined by the strong and electromagnetic interactions among the kaons, 
baryons and leptons. The time scales for these interactions are much smaller than those 
for weak reactions. Therefore, in setting up the kinetic equations, that determine the 
relatively slow change of chemical composition via weak reactions, the system may be 
assumed to be in thermal equilibrium, maintained by much faster strong and electromag- 
netic reactions. As a consequence, the physical quantities for the nonequilibrium state 
evolve adiabatically, adjusting to the gradual change of the chemical composition brought 
about by the weak processes. 

Assuming thermal equilibrium, we adopt the thermodynamic potential Q of the kaon- 
condensed phase which was derived by Tatsumi and Yasuhira from chiral symmetry ||38|| , 



including the thermal and quantum fluctuations around a condensate. In the present 
work, we shall neglect, for the sake of simplicity, both zero-point and thermal fluctuations 
in the thermodynamic variables except for the thermal contribution to the kaon number 
density, which allows us to see the effects of thermal kaons on the kinetics of condensation. 
Thus our expression for Q corresponds to the heavy-baryon limit for the nucleons ||38|| . 
With these assumptions, the thermodynamic potential per unit volume Q/V oi the kaon- 
condensed phase is written as 

n/v = {nc + ^K + nN + ne)/v , (6) 

where Qc and Qk are the potentials arising from a condensate and thermal fluctuations for 
the kaons, respectively, and Qn and fie are the corresponding potentials for the nucleons 
and the electrons, respectively. Here, Q, the number density Ui and the chemical potential 



We use the units in which h=c=kB=l throughout this paper. 
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fii {i = p,n,e ,K) are related to each other by = —d{Q/V)/dfii. Specifically, the 
classical contribution from the condensate, fi^, is given by 



nc/V = /^m^(l - cos 6) - -/VI sin^ Q , 



(7) 



where mx is the free kaon mass 
Ox, is given by 



494 MeV). In contrast, the thermal kaon contribution. 



-[w+(px)+/ix]/T 



(1 



— e 



with cc;+(p;^) [a;_(pii-)] the excitation energy for (cf. Sec. PB| for the explicit 

expressions for uj±{pk) ), and the nucleon contribution fi^v is written in the standard 
form: 



where the energy density for the nucleons is given by 
3 (37r2^2/3 



(9) 



N 



5 2m 



N 



+ nl/') - f{a + 26^,0(1 - cos^) + V,yUnB){np - ri„) /ub . (10) 



The first term in the r.h.s. of Eq. ([T0| ) is the nucleon Fermi-gas energy density, with 
rriN the nucleon mass, and the quantity a in the second term defined as cr = nsExN//^, 

and b = (up + -n„^/(2/^) is the y-spin density [see Sec. |IV A| for specific values]. The 

terms involving a and b, respectively, come from the s-wave K~N interaction given by 
the KN sigma term and the Tomozawa- Weinberg term. The expression for the energy 
contribution from kaon-kaon and kaon-nucleon interactions is essentially determined in a 
model-independent way by chiral symmetry. The third term in Eq.(]l0|) is the potential 
energy contribution to the symmetry energy. For V^ym(^B), we use the expression given 
by Prakash et al. p9|: 



5 J 



'111 



where 5*0 (= 30 MeV) is the empirical symmetry energy, and 6^,0 is the Fermi energy 
in the symmetric nuclear matter at the density no, with the function F{nB) taken as 
F{nB) = ub/uq for simplicity. The chemical potentials for the proton and the neutron, 
respectively, in are given by /ip = dSjy/dnp and /i^ = dSN/dun, respectively, so that 
the difference between fip and /i„ in the condensed state is written as 



iSn^npf/' (37r2n„)2/3 



2m 



N 



2m 



N 



+ 4:Vsym{nB) ■ {Up - nn)/nB 



-/ii^(l - COS I 



For f2e, we use the ultrarelativistic form for electrons. 



12) 



(13) 



47r2 127r2 ' 

where /le is the electron chemical potential and Ue [=fil/{3TT'^)] is the electron number 
density. Note that fie 7^ fJ^K when the system is not in chemical equilibrium. 
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B. Thermal kaon excitation 



As can be seen in Eq.(^, only the thermal kaon excitations constitute the thermal 
fluctuations in the thermodynamic potential fl, because they are responsible for the for- 
mation of a condensate through the weak reactions KT[(^)] . The excitation energy for 
the thermal kaons, uj±{pk), appearing in [Eq.® ], is obtained from the dispersion 
relations for kaon modes in the condensed state. The expression for uj±{px) depends on 
the method for treating the fluctuations around the condensate but numerical 

results have shown that there is little difference between the two methods used in Refs. 
5E| and 1^ . Later on we use the result of Ref . : 



u;±{pK)^±{b + fiK {cose- l)}+[p% + {b^ + mf) 



1/2 



(14) 



where m*^ = iri§ cos 6 = {fn\ — a) cos 6, and mj^ denotes the effective kaon mass which 
is reduced due to the KN scalar interaction Skn HI- Using (0) one obtains the total 
kaon number density from the relation, uk = —d{Q/V)/dnK'- 



riK = Ck + rij^ . 
Here (k is the condensed part of the kaon number density : 

CK = ~d{Qc/V + QN/V)/dfiK 
= fiKf sin^ + '2b f{l - cose) , 



(15) 



(16) 



and 71^ is the thermal part: 



— cose J d^PRfxiPK] 



with 



/i^(Px) 



i2ny 



1 



1 



(17) 



g[(^_(pK)-MK]/r _ I g[<^+(px)+MK]/T _ 1 



'1^ 



where the first and second terms are the Bose-Einstein distribution functions for the 
K~ and mesons, respectively. It is to be noted that the expression (|l^ is slightly 
different from the usual one for the noncondensed state by a reduction factor cos6' due 
to the existence of the condensate. The expression ([T6| ) is equivalent to the strangeness 
number density (k = {K~\S\K~) with the strangeness operator S = 2{Q — I3 — B). 
By the use of the classical field equation for e, dVt/de = PJ38[], 



sin e ( /i^ cos e + 2bfiK 



m 



*2 
K 







(19) 



which is valid where fluctuations of kaons can be neglected , it can be seen from ([l^) 
that, in the condensed phase {e ^ 0), the lowest excitation energy of K~ is equal to 
the kaon chemical potential, i.e., a;_(px = 0) = fix- This soft mode results from the 
spontaneous V^-spin symmetry breaking in the condensed phase (Goldstone mode), and 
leads to a divergence of fxiPK) at Pk=0 in the condensed state. On the other hand, 
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in the limit ^ = in Eg. ([M]) , one obtains the excitation energy of kaons in the normal 
(noncondensed) phase: 



±b 



m 



K ) 



1/2 



(20) 



In this case, there is a gap between the lowest excitation energy u;_(p/^ = 0) and ^k, so 
that fxiPK) has no singularity. The appearance of a soft kaon mode is also related to 
an instability of the normal state with respect to the onset of a condensate. (See also 



III. KINETICS OF CONDENSATION 
A. Rate Equations 

The thermally excited kaons, which are produced via the weak reactions KT, are con- 
verted into a condensate through kaon-nucleon and kaon- kaon scatterings: KN — >■ {K)N, 
KK K{K). The time required for the conversion of the thermal kaons to a condensate 
through strong interaction collisions is negligible compared with the time scale governed 
by the weak reactions. This is consistent with the assumption that the system is in 
thermal equilibrium, as discussed in Sec.[ll7\|. Thus, we put aside the detailed conversion 



mechanisms of the thermal kaons into a condensate, and study the kinetics of condensa- 
tion due to the weak reactions to obtain the typical time scales for the development of a 
condensate. With the kaon number density given by Eg. ([T5|) , one can see that the con- 
densate appears spontaneously once the number density of the thermal kaons is saturated 
for a given temperature. With this simplification, we can describe the growth of the con- 
densate during the whole nonequilibrium stage by following the kinetics of the condensate 
semiclassically, avoiding the discussion of quantum nucleation of the condensates. 

In general, the heat released by dissipation of energy in the system is expected to in- 
crease the temperature T. The temporal change of the temperature, T{t) can be obtained 
from the equation that gives the rate of change of the internal energy of the system: e.g., 

dE/dt = -e^ - ep (21) 

with e^, ep the luminosities of neutrinos and antineutrinos, respectively, for the neutrino- 
free-streaming case. For simplicity, however, we take the temperature to be constant 
during the nonequilibrium process. A more realistic calculation including the variation of 
the temperature will be discussed in the future. 

The number density of each chemical species is determined by the rate equations, 
with the rates of change of the electron (rig), the kaon {uk), and the proton (up) number 
densities given by 

+ r(^^U-F)(^(MU)^^)^y) _ r(MU-B)(^(MU)(^)^y) ^ (22a) 

+ r(KT-F)(^(KT)(^)^y) _ r(KT-B)(^(KT)(^)^y) ^ (22b) 
dnp{t)/dt = r(^^U-F)(^(MU)(^)^y) _ p(MU-B)(^(MU)(^)^2^) 
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(22c) 



where r'-"-' are the reaction rates («= KU, MU, KT), and 'F' ('B') denotes the for- 
ward (backward) process. As a consequence of baryon number conservation, the rate of 
change of the neutron number density is determined by dnp{t)/dt through the equation 
dnn{t)/dt = —dnp{t)/dt. 



B. Reaction rates 



Since the kaon-condensed state is obtained from the underlying chiral symmetry 
|37|, the matrix elements for the relevant reactions can be calculated from the chirally- 



transformed weak Hamiltonian. The expressions for these reaction rates have been given 
in paper I |^^. Here we only show the results. For the forward KT process (pa]), denoted 



as KT-F, one obtains 

r(KT-F)(^(KT)^y) = -J^Lcpf /sin^CCOS^^CCOS^^I 'P^^''^'' 




9(27r)7V^" ' ^ ^ 2|p^(n)|2 + m2 
X (m:)3^;T5/(KT)(^(KT)^y) ^23a) 



.,,0X10»)(^)(^) (|.)™4r|.-,,K..„ 

(cm-3 ■ s"^) , (23b) 

where = (/i^ + /^p — /^n) /T and I^^'^^u^ T) is the integral over the kaon momentum 
divided by T, with x = \pk\/T: 

^' ^ 6 Jo (cul(x) +/ii^/T)3 1 - e-^-W e^-{x)+u_i ^ > 

where ul{x) = {uj_{x)— hk) IT . In ( p3a| ), 5'^(=1.25) is the axial- vector coupling strength, 
/ = /vrATAr/m^ the IT NN coupling strength divided by the pion mass, Gp the Fermi cou- 

Q 

pling constant, and ^c'(~0.24) the Cabibbo angle. The factor 5f^G'i?/sin^c'Cos6'ccos^ - 
originates from a npK~ vertex factor in the transition matrix elements obtained from the 

Ipf('"')P 

chirally-rotated hadron currents, whereas the other factor, j- with prii) 

Ipf(^)| + 'i^i 

{i = n, p) the Fermi momentum, comes from the one-pion exchange potential introduced 
to describe the long-range part of the interactions. Due to the approximate kinematical 
condition ||^ , only the nucleon axial- vector current (proportional to qa) is found to con- 
tribute to the KT reactions after all the matrix elements for the lowest-order diagrams 
are summed. Hence the current for the s-wave condensate, which has only a time com- 
ponent, does not couple to the nonrelativistic nucleon current. The remaining factor in 
( p3a|) arises from the phase space integrals in which m|^f is the effective nucleon mass and 
Tg, the temperature in units of 10^ K. For simplicity, we take the values of the nucleon 
effective masses to be m*/mAr = m*/mAr=0.8 |]32| . 

For the forward KU process (KU-F) [(|3a|)], one obtains 



r(KU-F)(^(KU)^y) = -g^sin2^csin2^{lO + 3(^i + 9^1)HVeT^/2(e(^^) 
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= (6.6 X 10^9) ( !!!kyj^ sin2 ^(t)T95/2(^(KU)) ^^^-3 . ^-i^ ^25) 
where gA = F -\D = 0.15 with F + D = = 1.25, D/{D + F) = 0.658 |2g, hiu) = 

roo 

f dxx'^[Tr'^ + (x + 'u)^]/(l + exp(x + u)) and ^(^u) = (/ig — fij^)/T. That the rate for 
Jo 

KU-F is proportional to sin^ 9 is derived from the fact that the matrix element for KU-F 
is proportional to sin^. As a consequence, the KU reactions are operative only when a 
condensate is present. 



For the forward MU process (MU-F) [(^)], we refer to Haensel's result [41 



p(MU-F)(^(MU)^y) = (5.9 X 1023)f!^y^'cos2^Tjj2(e(''^)) (cm-^ ■ s-^) , (26) 

POO 

where J2(m) = / rfxx^[97r^ + 10tt'^{x + + (x + m)V(1 + exp(x + u)), and ^(^u) ^ 
Jo 

{Hp + He — Hn)/T. Here, the matrix elements have been slightly modified by inclusion of 
an additional factor of cos^(6'/2) due to the presence of a condensate (cf. Paper I). 

The backward processes, denoted by the suffix B, are related to the forward processes 
by way of the following relations: 

r(KT-B) (^(KT) ^ y ^ _ ^^(KT) ^(KT-F) ^^(KT) ^ ^ (27) 

and 

p(KU-B)j^^(KU) 2^-) ^ p(KU-F)j^_^(KU) j^-j pCMU-B) j^^(MU) j^-j ^ -[^(MU-F) j^_^(MU) rp^ (28) 

which are valid at low temperatures The relation ( ^7] ) between the forward and 

backward reaction rates for KT (where bosons are involved) differs from the corresponding 
relations (^Sf ) in the case of the KU or MU process (where only fermions are involved). 

C. Initial conditions 

We adopt the following initial conditions: At t=0, we assume normal neutron-star mat- 
ter (6 = 0), composed of nonrelativistic protons (p), neutrons (n), and ultra- relativistic 
free electrons (e~), which is charge neutral, = n^. The baryon number density Ub 
for the noncondensed state is taken to be larger than the critical density for kaon con- 
densation, (=3—4 no). Initially, the system is assumed to be in /^-equilibrium, i.e., 
/j,^ = Hp + He [ ^^^^\t = 0) = 0] due to the rapid /3-decay reactions, n ^ p + e~ + Ue, 
p + e~ ^ n + Ue, when the neutron star is hot at an early stage. The initial values for the 
proton fraction Xp = Up/ns and for the electron chemical potential He then obtained 
from the charge neutrality condition, ubx^ = (/ig)^/(37r^), and the /^-equilibrium condi- 
tion where Eq.(|l^ is used with 9 = 0. Finally, the total strangeness is assumed to be 
almost zero corresponding to equal numbers of thermal iT+'s and ^^'''s present, which, 
using Eqs. (0), ( 1T7|) and (|^) with 9 = 0, gives the initial value of the kaon chemical 



^The superscript '0' denotes the initial value at t = 0. 
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potential /z^c = —iP- As a consequence of these conditions, the initial kaon chemical po- 
tential /i^f has a large negative value, while /ig is positive. Thus the system is far from 
chemical equilibrium with |^(KT)(t = 0)| = \i^^^\t = 0)| = - iJ,%)/T\ > 1. For 
example, one obtains ^^^^\t = 0) = -46 and C^^^\t = 0)=46 for J:kn= 300 MeV, nB= 
0.55 fm~^ and T = 1 x lO^^K. Starting from these initial conditions, one can first obtain 
the number densities ni{t) at later times from the rate equations (|2^). Next one can 
obtain the electron chemical potential /ie(^) from the relation yUe(^) = (37r^r7,e(t))^^'^, the 
kaon chemical potential fixit) and the chiral angle 9{t) from Eqs.(p!5D — (p!8D and (|1^), and 
then the difference between ^p{t) and finit) from Eq.(|T^). Because the charge neutrality 
condition, Upit) = rixit) + rieit), is built into the rate equations (|2^), it can be seen that 
this condition is always satisfied. The system evolves dynamically toward an equilibrated 
kaon-condensed phase through the nonequilibrium weak processes, KU, MU, and KT. 



IV. NUMERICAL RESULTS AND DISCUSSION 

We choose the value for the KN sigma term to be S/^Ar=300 MeV The critical 



density n§ is then estimated to be n§=0.49 fm^^ (= 3.0 uq ) [c.f. Ref. ^3|]]. In Table 
|, we list the values for /i^^ and the proton fraction Xp for the initial noncondensed state 
(t=0), and those of 6, fix and Xp for the kaon-condensed state in chemical equilibrium 
{t — > oo) [which we denote by the superscript 'eq']. It has been shown in Ref. that 
the temperature dependence of the critical density and physical quantities in the kaon- 
condensed phase is weak for the temperature less than several tens of MeV. Thus all the 
values listed in Table | have been estimated at T=0. Here, we study the temporal behavior 
of the physical quantities for two different baryon number densities, 71^=0.55 fm~^ and 
0.70 fm~^. The density 72^=0.55 fm~'^ is close to such that we have a rather weak 
condensed state with small 6'^^. For nB=0.70 fm~^, on the other hand, we have a well- 
developed condensed state with a large order parameter O^"^. When the baryon number 
density is increased above the critical density for kaon condensation, the number density 
of negatively- charged kaons increases as the condensate develops. Charge neutrality of 
the system is maintained by the increase in the proton number density. Consequently, the 
proton fraction Xp^ increases as the baryon number density increases in an equilibrated 
kaon-condensed state. On the other hand, the negatively charged electrons are replaced 
by the negatively charged kaons as a condensate develops with an increase in the baryon 
number density. The consequent decrease in the electron abundance produces a decrease 
in the charge chemical potential fif^ (=/i^), resulting in a large proton fraction and a 
reduced charge chemical potential, which are two characteristic features of the kaon- 
condensed state 



A. Two typical time scales 

In Fig.l (a), we show the chiral angle 5 as a function of time for 72^=0.55 fm" 



-3 



and T=1.0xl0^i K. For the initial values, /x?^ = -6° = -139 MeV and a;2=0.14. The 



p 



evolution of the kaon condensate may be divided into the following three stages: (I) 
no condensate {0{t) = 0) with thermal kaons present until the onset of condensation, 
(II) onset of condensation and its monotonic growth, and (III) the asymptotic stage 
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near chemical equilibrium. In equilibrium, one finds ^^^'^=0.48 (rad) from Fig. |T|(a), and 
the kaon chemical potential and the proton fraction have the values /i^=203 MeV and 
= rf^/n-Q = 0.23, respectively (see Table 1). 

The onset of condensation from a noncondensed state and its subsequent growth are 
related to the change in shape of the thermodynamic potential as a function of 6 around 
6 = 0. This change in shape is also associated with the appearance of a soft kaon mode, 
reflected in a singularity of the distribution function fxiPK), Eq-(0)- To see this, we 
expand fl{t)/V around ^ = 0: 

m/v = -^D],'eitr+o{eitr) , (29) 

with D]^^ = D]^^{u! = ^K(t),pK = O-jUb), where D]^^{u,pK',nB) is the inverse kaon 
propagator: 

Dj,\u, pk, ub) =u^-pI + 2b{t)uj - . (30) 

Beginning with = —b^ at t=0, the kaon chemical potential fJ-xit) increases mono- 
tonically with time. As a result, the factor D]^^ changes sign, depending on the value 

1 /2 

of hk relative to the critical value, /i^(t) = -b{t) + {b{tf + mf) , defined as a root 

of D]^ . In Fig. 1^, we show the thermodynamic potential per particle Vt/iVuB) as a 
function of 9 for nB=0.55 fm~^ and T =1.0x10^^ K and the other parameters and /ij 
{i = p,n,e~ , K~) fixed at three times corresponding to cases (i) — (iii). Case (i) corre- 
sponds to a time before the onset of condensation, for which fi^ is smaller than /i^. [E.g., 
at t=0, Hk = —b^ < fJ'xi^ — 0)-] this case, D]^^ is negative, which means that the ther- 
modynamic potential per unit volume Q/V is convex at ^=0, as seen from Eq. (|29| ) and 
the curve (i) in Fig. ^, and the noncondensed state (6'=0) is stable against fluctuations 
in 9. Because /i^ corresponds to the lowest excitation energy of the K~ , U-{pk = 0) m 
normal matter, one can see that uj^{pk) > uj^{pk = 0) = yU^^ > fiK- Thus the Bose- 
Einstein distribution function fx here has no singularity. Case (ii) corresponds to a time 
for which fixit) = fixit): where D]^^ = 0. In this case, the system becomes unstable with 
respect to the fluctuations in 9 (the curve (ii) in Fig. ||), and since u_{pk = 0)(= /i^) 
is equal to fix, fKiPx) becomes divergent at Pa'=0. Thus the appearance of a soft kaon 
mode with uj^{pk = 0)=fiK implies the onset of an instability due to the formation of a 
condensate. Finally, case (iii) corresponds to a time after the onset of condensation, for 
which yUi^(t)exceeds /U^(t), so that D]^^ > 0. In this case, fl/V is concave near 6^=0 as 
a function of 9; which means that the noncondensed state corresponds to the maximum 
of ^l/V, as seen from the curve (iii) in Fig. |^, and is unstable against the formation of 
a condensate. The system therefore evolves into a condensate described by 9{t), which 
is determined adiabatically by a minimum in the thermodynamic potential [cf. Eq.(p^)]. 
By the use of the dispersion relations for the kaon in the condensed phase (see |1I B]) ||38| , 
it can be shown that the lowest excitation energy of in the condensed phase is equal 
to fiK, i.e., lj_{pk = 0) = fiK(t); there is a soft kaon mode in the condensed phase, 
and fxiPK) is divergent at pk = 0, as in the case (ii). From Fig. 1(b), where D]^ is 
shown as a function of time, one sees that the function D]^ changes sign at the onset of 
condensation. We can call the onset time for condensation the nucleation time, Tnuci, and 
the time for growth of the condensate the coherence time, Tcoh, after Stoof's discussion 
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about the formation of BEC in an atomic gas ^M. From Fig. P, one finds Xnuci = 8 x 10 



sec, and Tcoh = 5 x 10 ^ sec for nB=0.55 fm and T = 1.0 x 10 K. 



B. Connection between the reaction rates and the change in chemical composition 

The temporal behavior of the reaction rate for each weak process is shown for nB=0.55 
fm~^ and T = 1.0 x 10^^ K in Fig. ^. The sohd, dashed and dotted hnes denote the KT, 
MU, and KU processes, respectively. The curves show that the forward thermal kaon 
process KT-F is dominant in magnitude over the KU and MU reactions throughout the 



equilibration process [0. Hence the production of the thermal and condensed kaons 
proceeds mainly via the KT-F reaction, which is responsible for the onset of condensation 
and its subsequent growth. 

In Fig. ^, we show the temporal behavior of the dimensionless parameters ^'-"^ {a= KU, 
MU, KT), which measure the deviation from chemical equilibrium. Beginning with a large 
negative value, increases monotonically toward zero, where the system is in chemical 
equilibrium with respect to the KT reactions. Note that the KT reaction rates are large 
even in the initial stage where the system is far from chemical equilibrium ^ 1)- 

In Paper I, we examined in detail the characteristic roles of the soft and hard thermal 
kaons in the KT reactions, and demonstrated that the high-energy component of the 
thermal kaons contributes to the KT reaction rates during the noncondensed stage. On 
the other hand, a low energy component (a soft mode ) of the thermal kaon excitations in 
the kaon-condensed state was shown to contribute mainly near chemical equilibrium ||37 |. 



Let us now examine how these reactions change the chemical composition. The number 
densities of the chemical species are shown in Fig. ^ for the same density and temperature 
as in Figs. |1]— In stage I (0 < t < Tnuci), where there is no condensate, the kaons 
produced by the KT reactions are thermal, and the KU reactions cannot proceed because 
there is no condensate {6 = 0). The initial /3-equilibrium is disturbed by the proton 
excess, created through the KT-F reaction, and the MU reactions are enhanced because 
the MU-B reaction proceeds more rapidly than the MU-F reaction, which restores /3- 
equilibrium by reducing the proton excess. Nevertheless, the proton number density 
still increases, since the KT-F reaction rate is larger than the MU-B reaction rate, as 
can be seen in Fig. ^ (see also Paper I). Because the total baryon number density n-s is 
constant, the neutron number density n„ therefore decreases with time. In stage I, the 
change in the electron number density rif, is brought about only by the MU reactions, 
and rif, decreases slightly with time, the MU-B reaction proceeding more rapidly than the 
MU-F reaction. The extent of the changes in Hp, n„, and Ug is very small. Consequently, 
although the electron chemical potential fie decreases with time because of the decrease 
in He, this potential, along with the proton and the neutron chemical potentials fip and 
fin remain almost constant in stage I. On the other hand, the kaon chemical potential fix 
increases significantly with time; fj}^ = —IP at t=0 and — jJ^^K at t = t^ucI- Hence 
which is proportional to /ie — f^K, decreases monotonically with time and increases 
monotonically with time [cf. Fig. while ^(^u) j-gniains close to its initial value (~ 0) 
except for the time near Tnuci- 

In stage II (rnuci <t< Tcoh), the KU reactions can proceed due to the appearance of a 
condensate, but the qualitative behavior of each chemical composition is similar to stage 
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I. In particular, the kaon and proton number densities increase through the KT reactions, 
while those of the electrons and neutrons decrease. The change in the number density 
of the kaons uk comes mainly from the condensed part (k [(IH)] after the saturation of 
the thermal part n^, with the condensed part (k roughly proportional to the square of 
6, (k oc 6'^, for a small chiral angle, as seen from Eq. ([T6|). According to Fig.[^(a), the 
chiral angle 6 increases rapidly with time for t ^ Tnuci- Thus the change in the number 
density of the kaons rix is substantial, being more pronounced in stage II than in stage I 
due to the appearance of a condensate. The change in the number density of the electrons 
He in stage II is caused by the KU and MU reactions. But since the reaction rates for 
the electron absorption processes, KU-B and MU-B, are smaller by orders of magnitude 
than the rate for the kaon production process KT-F, the change in the number density 
of the electrons is less marked than the change in the number density of the kaons, and 
remains almost constant near the onset of condensation, until both the MU-B and KU-B 
reactions become maximum and operate significantly to reduce ng in the later part of 
stage II (cf. Figs. ^ and From the temporal behaviors of uk, n^, the charge neutrality, 
np = He + rtKi and baryon number conservation, rip + ?t,„ = n^, it can be seen that the 
change in the number densities of the protons and neutrons become remarkable in stage 
II. 

During the earlier part of stage II where the number density of the electrons is almost 
constant, the parameter ^(^^^ remains unchanged because both the electron and the kaon 
chemical potentials are almost constant. On the other hand, as a result of the increase in 
the proton number density and the decrease in the neutron number density, the difference 
between the proton and neutron chemical potentials becomes larger with time, as seen 
from Eq.([T^). Thus the parameter ^(^u) increases with time, after which both and 
^(MU) (igCTgasg -with time, and the system enters into the final stage III (cf. Fig. p. 

The coherence time Tcoh may be identified with the time at which the magnitudes of 
the forward and backward reaction rates for KT become equal. As seen in Fig.^, chemical 
equilibrium with respect to the KU and MU reactions is achieved at a later time than 
that at which the KT reaction reaches equilibrium. 

Finally, in stage III [t > Tcoh), where the system is close to chemical equilibrium 
^^(KT) ^ g^^ both the KU and MU reactions compete with each other to determine the 
dynamical behavior of the system. In this stage, the deviation parameters ^(^^^ and 
^(MU) fianip exponentially with time, and the system approaches chemical equilibrium. 
The relaxation time Xrei can be calculated analytically [see Eq. ( |A11| ) in Appendix A]. 
For ^3=0.55 fm~^ and T=1.0xlO^^ K, the analytic result gives Trei = 2 x 10~^ sec, which 
is in agreement with the numerical result read from Fig. ^ 



C. Temperature-dependence of the characteristic time scales 



Next we compare the temperature dependence of the characteristic time scales, Tnuci, 
Tcoh, and Trei at a fixed baryon number density. Figure || shows the behavior of the chiral 
angle 9 for nB=0.55 fm~^ at three different temperatures, T = 1.0 x 10^*^ K, 1.0 x 10^^ 
K, and 5.0 x 10^^ K. In addition. Table ^ lists the values of Tnuci and Tcoh, estimated 
from the numerical calculations, and the value of Tfci, obtained from Eq. ( |A11|) . Roughly 



speaking, it can be seen that Tnuci depends weakly on temperature, while Tcoh and Tj-ei 
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depend sensitively on temperature: In particular, from Table |T| and Fig. one obtains 
Tnuci=2 X 10~^ sec 6 X 10"^ sec, rcoh=2 x 10~^ sec ^ 4 x 10"^ sec, and rrei=40 sec 
5 X 10"^ sec as T = 1.0 x 10^° K ^ 5.0 x 10" K for nB=0.55 fm'^. 

Consider first the temperature dependence of the nucleation time, Tnuci, which provides 
a typical time scale for saturation of the thermal part of the strangeness number density 
n^(t). For t > Tnuci, n^i^) changes little. In chemical equilibrium, the thermal kaons 
occupy a progressively larger part of the total strangeness density as the temperature 
increases, and two competing effects produce the temperature dependence of Tnuci- (i) 
The value of n^[t) at saturation is larger at higher temperature, which tends to increase 
the nucleation time, Tnuci- On the other hand, (ii) a higher temperature gives a higher 
reaction rate for the relevant KT process, (cf. Paper I). As a consequence the thermal 
kaon number density n^{t) comes to saturation earlier. 

In contrast, the temperature dependence of the coherence time Tcoh is explained by 
the fact that after a condensate appears, its subsequent growth in stage II is controlled 
mainly by the KT reactions, as seen in Sec. |IV B| , and the KT reaction rates depend 
sensitively on the temperature over the entire nonequilibrium processes, (cf. Paper I). 
This leads to a sensitive temperature dependence of Tcoh- 

Finally, the relaxation time Trei depends on the KU and MU reaction rates through 
the quantities f^^^-F) ys^, p(MU-F) ^7^, ^^^^ ^q. in Appendix A], and 

from Eqs. ( [All| ) and ([A9|) , it can be seen that the relaxation time Trei depends sensitively 
on the temperature, so as to decrease significantly as the temperature increases. 



D. Density-dependence of the characteristic time scales 

Here we discuss the temporal behavior of the system at a higher baryon number 
density, nB=0.70 fm~'^, for which one can expect a fully developed kaon-condensed phase 
when the system is in chemical equilibrium. Initially we start with nj^{t = 0) = 0, which 
gives /i^ = —W = —180 MeV and Xp=0.16. The temporal behavior of the chiral angle 6 
is shown in Fig. |^ for nB=0.70 fm~^ for the same three temperatures as in Fig. The 
numerical values for the nucleation time Tnuci, the coherence time Tcoh, and the relaxation 
time Trei for the density and temperatures corresponding to the three curves shown are 
listed in Table fT[ From an examination of this table and Fig. 0, it can be seen that while 
the nucleation time Tnuci depends weakly on the temperature, the coherence time Tcoh and 
the relaxation time Tjci depend sensitively on the temperature, becoming shorter as the 
temperature is raised. These features for Tuuci, Tcoh, and Trei are qualitatively the same as 
in the case of the lower baryon number density nB= 0.55 fm~^. 

The time scales obtained for nB=0.70 fm~^ can be compared with those for nB=0.55 
fm~^. From Table |l[ one can see that the nucleation time Tnuci is sensitive to the baryon 
number density, with the ratio Tnuci (0.55 fm~^)/Tnuci (0.70 fm~^) equal to 2x10 — 2 x 10^ 
for 1 X 10^° K< T < 5 X 10^^ K. This density dependence is explained as follows: In 
the noncondensed state where l^^^"*"^] S> 1 (^*^^"'") < 0), the KT-F reaction rate is very 
sensitive to l^^^"""-*] and becomes large rapidly as |.^*-^"'"'*| increases [0. In addition, the 
kaon chemical potential fix (<0 at t = 0) in |^(^'^)| increases rapidly with time in stage I, 
and, before the onset of condensation, it is already almost equal to the equilibrium value 
A^ic (>0), which has a smaller value for a larger baryon number density. As compared 
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with the density dependence of /j^k, the density dependence of the difference between the 
proton and neutron chemical potentials fXp — /i^ (< 0) in |^*-^"'"''| is weak. Hence, as the 
baryon number density increases, the parameter l^*-^"""^! (=|/ii^ +/ip — /i„|/T ^ 1) increases 
mainly due to the decrease in ^x- As a result, the KT-F reaction rate in stage I is much 
larger for a larger baryon number density, which considerablly reduces the time required 
for saturation of the thermal kaon number density. 

After the appearance of a condensate, the value of I^C^^"'"''! becomes sufficiently small, 
so that the density dependence of |^(^'^)| becomes insignificant. Therefore, the KT-F 
reaction rate does not depend much on n^. Hence the coherence time T^oh, which is 
controlled by the KT-F reaction, shows only a weak dependence on ub- In a similar way, 
near chemical equilibrium where |^*-^^''| <^ 1 and l^*-^^-*] -C 1, the density dependence of 
both the KU-F and MU-F reaction rates is weak, and therefore Trei, which is determined 
by these reaction rates, is found to have a weak dependence on n^. 

E. Implications for the delayed collapse of a neutron star 

From Table |I|, one can see that a condensate develops fully with a characteristic 
time scale given by the coherence time Tcoh- In the context of the delayed collapse of a 
newly born neutron star, we consider some implications of the results. Since the neutrino 
degeneracy is not taken into account in our framework, our results are applicable to the 
cooling era after the deleptonization. Hence the delayed collapse associated with a kaon 
condensate should be considered in this cooling era. 

The low temperature case (e.g., T = 1 x 10^° K) may apply to the final stage of the 
cooling era, while the high temperature case (e.g.,T ^ 10^^ K) to the initial stage. In 
both cases, we have found that the time scale of Tcoh is very small; Tcoh ~ 0.1 sec for 
T = 1 X 10^° K and Tcoh ~ 5 x 10~^ sec for T ^ 10^^ K, which should be compared 
with the cooling time scale of order of ten seconds. Hence we can conclude that, in both 
cases, the time delay of a collapse due to the formation of a condensate is negligible as 
compared with the cooling time scale and will have a minor effect on the evolution of the 
protoneutron star. 

V. SUMMARY AND CONCLUDING REMARKS 

We have considered the kinetics of kaon condensation by the use of rate equations 
which include the three weak reactions: the thermal kaon process, the kaon-induced 
Urea process, and the modified Urea process. The thermal kaon process is shown to be 
dominant over other weak reactions throughout the equilibration process. The evolution 
of the kaon condensate is divided into the following three stages: (I) the noncondensed 
stage before the onset of condensation, (II) onset of condensation and its growth until its 
saturation, and (HI) the asymptotic relaxation stage near chemical equilibrium. The role 
of thermal kaons and especially the connection between the existence of a soft kaon mode 
and the instability of the noncondensed state has been clarified. It has been found that a 
full development of a condensate is characterized by a time scale given by the coherence 
time Tcoh- Using these results, we have made a brief comment on implications of our 
results on the delayed collapse mechanism of a newly-born neutron star into a black hole. 
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The situation adopted in our framework may be applicable to the cooling era after the 
deleptonization. It is found that the time scale of the development of a condensate (the 
coherence time) is much smaller than that of the cooling for the relevant temperatures 
(1 MeV— several tens of MeV). Therefore, the time delay for the formation of a kaon 
condensate may not affect the delayed collapse of a neutron star. 

Several effects which should be taken into account within our framework in the future 
are as follows: 

(1) We have used the thermodynamic potential neglecting any fluctuations except for those 
produced by thermal kaons. However, for a neutron star just born in a supernova, since the 
initial temperature is as high as several tens of MeV, thermal fluctuations other than the 
kaon thermal loop contributions |^ should be fully incorporated into the thermodynamic 
potential. 

(2) Furthermore, since the matter is opaque to the neutrinos at high temperatures, the 
thermal effects associated with neutrino diffusion may also be important. In particular, 
the degenerate neutrinos will contribute to the thermodynamic potential and influence the 
relevant reaction rates. Therefore, the neutrino diffusion may affect the thermal evolution 
of a newly-born neutron star. 

(3) The reaction rates for the relevant weak processes have been obtained in the low 
temperature approximation (c.f. Paper I). The expressions for the reaction rates must be 
extended to the high temperature case, and the phase-space integrals should be performed 



numerically beyond the low temperature approximation [46 



(4) During the equilibration process, we have assumed that the temperature T is constant: 
the energy produced by the weak reactions is lost to the surroundings. However, because 
the energy release due to the nonequilibrium weak reactions will heat up the matter, the 
temperature is in general time-dependent during the course of the equihbration process, 
and this can affect the evolution of the star. This thermal effect can be included through an 
additional equation that determines the rate of change of the internal energy of the system. 
Q A signiflcant increase in the temperature could possibly lead to further reduction of 
the time scale for the growth of a condensate represented by the coherence time. On 
the other hand, as a result of the raise in the temperature, the critical density for 
kaon condensation can be pushed up to higher density, which may cause the delay in the 
onset of condensation. In this respect, it is necessary to obtain a phase diagram in the 



T — ub plane ||38| , ^ , and the dynamics of the onset and growth of condensation must be 



considered along a trajectory in this plane. 
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APPENDIX A: ASYMPTOTIC BEHAVIOR NEAR CHEMICAL 

EQUILIBRIUM 

We derive an expression for the relaxation time which characterizes the asymptotic 
behavior of the system near chemical equilibrium (stage III). 

1. Linearization of ^(^^^t) and C^^^HO 

In stage III, both the forward and backward KT reaction rates are almost equal, i.e., 
^(KT)^^^ ~ 0. Therefore, we consider the case where the relaxation proceeds only through 
the KU and MU reactions, neglecting the KT reaction. The parameters $,^^^\t) and 
^^^^\t) are written as 

^^'^''\t) = [6fi,{t)-6f,K{t)]/T , (Ala) 
e("'U)(t) = [5/ie(t) + %(t) - 5/i„(t)]/T , (Alb) 

where 6ni{t) [i = p,n,e~ , K~) is the deviation of the chemical potential from the equi- 
librium value, i.e., 5fii{t) = jj,i{t) — nl'^. From the relation, /ig = (STr^ng)^/^, one obtains 

5/ie(t) = (7r//i:<^)2<5ne(t) , (A2) 

where Sne{t) = ne{t) — nl^. One can also write the deviations of the chemical potentials 
SfiK^t) and Sfipit) — Sfinit) in terms of 5ni{t) {=ni{t) — rf^) from linearized forms of the 
three equations with respect to 6fiKit), ^^it), Sfip(t) — 6finit), and Sni{t) represented by: 
the expression for the kaon number density [Eq. (0)], the classical field equation [Eq. 



(|19|) ], and the expression for the difference between /Xp and /i„ [Eq. (p^ j . Note that 
we can neglect the deviation of the thermal kaon part in the kaon number density, i.e., 
driKit) = SCxit) ((5n^(t) = 0), because the thermal kaon number density n^^t) is almost 
constant after a condensate appears until the system reaches chemical equilibrium. 

The deviation of the electron number density Sne{t) is divided into two contributions: 

6ne{t) = 5(^")ne(t) + 5(^^)ne(t) , (A3) 

where 5^^^^ne{t) [5^^^^ne{t)] is caused by the KU reactions (the MU reactions). The 
number densities for the other species are related to 5^^^^nf,{t) and 5^^^^ni.{t) by 

SuKit) = -5(^^)ne(t) , (A4a) 
Sripit) = ^ (A4b) 

6nn{t) = -6np{t) = -6^^^^ne{t) (A4c) 
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Following the above results, the deviation parameters 
of 6^^^^ne{t) and 6^^^^ 71^(1). The result is 



P 



T]) can be expressed in terms 

(A5a) 
(A5b) 



where 



P 



eq 
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8l^sym(y^B) 



eq 



In Eg . ( |A6|) , the first term comes from the deviation of the electron chemical potential 
5/ie, and the second term from the deviation of the parameters for a condensate, 

i.e., 6fiK and 66. The third and fourth terms in Eq. ([A6c| ) come from the deviation of 
the proton and neutron chemical potentials. It is to be noted that a condensate already 
develops fully in stage III, so that the quantities for the condensate are almost constant. 
Thus it is a good approximation to set 6fiK = 66 = 0, which leads to p = g 

3 2mj, )+ 



and r 



(A6a) 
(A6b) 
(A6c) 



2. Differential equations for ^(^^\t) and ^^^^\t) 



By differentiating both sides of Eq. ( |A5|) with time, one finds that the rates of change 
of ^(KU)^^^ and C^^^\t) are related to the reaction rates for KU and MU, through the rate 
equations 

l5(KU)^^(i) = r(KU-F)(^(KU)(^)^^) _ r(KU-F)(_^(KU)(^)^^) ^ (A7a) 
i^6^^^)n,{t) = r(MU-F)(^(MU)(^)^y) _ r(^U-F)^_^(MU)(^)^y) _ (A7b) 

For the system close to chemical equilibrium, ^ 1 and I^C'^^^H^)! 1- ^^is 

case, the r.h.s. of Eqs. (|A7a|) and ( |A7b|) can be expanded up to first order in ^^^^\t) and 
^(MU)^^^^ respectively. Thus one obtains the coupled differential equations for ^^^^\t) and 

A^(KU)(i) = _|f (KU-F) . ^(KU)(^) _ If (MU-F) . ^(MU) ^ (^Sa) 
l^(MU)(i) = _|f (KU-F) . ^(KU)(^) _ (MU-F) . ^(MU) ^ (A8b) 
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where 

p(KU-F) ^ 2ipMr(KU-F)(^(KU) = o,T) ^ 1.3 X r(KU-F) ^^(Ku) = Q, T) oc T^ (A9a) 
-'2(0) 



f(MU-F) ^ 2H^r(^u-F)^^(MU) = 0^2^) ^ 1.2 X r(^^u-F)(^(MU) = Q, T) oc T^ (A9b) 
-'2(0) 

with/2(0) = |7r2C(3) + f C(5)=41.13,/^(0) = (dhiu) / du) ^=0 = J2(0) = f A(3) + 

2257r2C(5) + ^C(7)=4598, and J^(0) = (dJ2(M)/ciu)„=o = -ffTi-^- From Eq. one 
obtains the following forms of the solutions: 

^(KU)(^) = Cie"*/"" + Cse"*/^ , (AlOa) 
^(MU)(^) = Cae-*/"" + C4e-*/"'' , (AlOb) 

where Ci — C4 are constants, and and r;, are time constants which characterize the 
asymptotic evolution of the system. The relaxation time T^ei is defined as the larger of 
and Tb, and one obtains 



ypf (KU-F) _^ ^p(MU-F) _^ r{pf(™-^) - rf(^U-F)|2 _^ 4^2p(KU-F)p(MU-F)] ^2 



Trel 



(pr - g2^r(KU-F)r(MU-F) 

(All) 
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TABLES 



TABLE I. Quantities for the initial noncondensed state {t = 0) and for the chemically equilibrated 
state {t x). The former (the latter) quantities are denoted by the superscript '0' ('eq'). All the values 

are estimated at T=0. 



nB(fm 


( MeV ) 


™0 


e^'i (rad) 


/^k' ( MeV ) 




0.55 


-139 


0.14 


0.48 


203 


0.23 


0.70 


-180 


0.16 


0.91 


114 


0.39 



TABLE II. The nucleation time Tnuci, the coherence time Tcoh, and the relaxation time Trei for 
different densities and temperatures. 



nB(fni ^) 


T(K) 


Tnucl (sec) 


Tcoh (sec) 


Trel (sec) 


0.55 


1x10^° 


2x10"^ 


2x10"^ 


4x10 




1x10^1 


8x10"^ 


5x10-5 


2x 10-4 




5x10" 


6x10-9 


4x10-^ 


5x 10-^ 


0.70 


1x10^0 


1x10-" 


6x10-2 


9x10 




1x10" 


3x10"^° 


1x10-5 


2x lO-'' 




5x10" 


3x10"^° 


2x10-^ 


Ix lO-'^ 
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FIGURES 




0.0 

10-' 
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10-7 10-6 10-5 10-4 

t ( sec ) 'f ^^'^ ■f 

'CnucI ^coh ^nucl ^coh 

FIG. 1. (a) Chiral angle 9 and (b) D^^, the negative of the expansion coefficient for the 
thermodynamic potential per unit volume i}/V with respect to the squared classical kaon field, 
as a function of time for the baryon number density nB=0.55 fm~^ and the temperature T 
=1.0x10^^ K. 



> 

0) 



III 

> 




1.0 1.5 2.0 

d ( rad ) 

FIG. 2. Thermodynamic potential per particle Q/ {Vns) as a function of 9 corresponding to the three 
cases (i)— (iii). The long-dashed and short-dashed lines correspond to the cases at t = and t ^ oo, 
respectively. The dots on the lines represent the locations of the minima of O. 
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FIG. 5. Temporal behavior of the number densities of the chemical species. 
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FIG. 7. Temporal behavior of the chiral angle 6 for nB=0.70 fm ^ and three different temperat 
= 1.0 X 10^° K, 1.0 X 10^1 K, and 5.0 x 10^^ K. 
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